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ABSTRACT 

Currently-proposed galaxy quenching mechanisms predict very different behaviours during 
major halo mergers, ranging from signihcant quenching enhancement (e.g., clump-induced 
gravitational heating models) to signihcant star formation enhancement (e.g., gas starvation 
models). To test real galaxies’ behaviour, we present an observational galaxy pair method for 
selecting galaxies whose host haloes are preferentially undergoing major mergers. Applying 
the method to central L* (IO'^'M© < M* < lO'^'-^M©) galaxies in the Sloan Digital Sky Survey 
(SDSS) at z < 0.06, we hnd that major halo mergers can at most modestly reduce the star¬ 
forming fraction, from 59% to 47%. Consistent with past research, however, mergers accom¬ 
pany enhanced specihc star formation rates for star-forming L* centrals: ~10% when a paired 
galaxy is within 200 kpc (approximately the host halo’s virial radius), climbing to ^ 70% 
when a paired galaxy is within 30 kpc. No evidence is seen for even extremely close pairs 
(< 30 kpc separation) rejuvenating star formation in quenched galaxies. For galaxy forma¬ 
tion models, our results suggest: (1) quenching in L* galaxies likely begins due to decoupling 
of the galaxy from existing hot and cold gas reservoirs, rather than a lack of available gas 
or gravitational heating from infalling clumps, (2) state-of-the-art semi-analytic models cur¬ 
rently over-predict the effect of major halo mergers on quenching, and (3) major halo mergers 
can trigger enhanced star formation in non-quenched central galaxies. 

Key words: galaxies: haloes; galaxies: formation 


1 INTRODUCTION 

In the context of Lambda Cold Dark Matter (ACDM) cosmologies, 
dark matter halo (gravitationally self-bound structures) masses cor¬ 
relate strongly with the stellar masses of the galaxies at their cen¬ 
tres (e.g.. More et al. 2009; Yang et al. 2009; Leauthaud et al. 2012; 
Reddick et al. 2013; Watson & Conroy 2013; Tinker et al. 2013). 
Indeed, one-to-one matching of galaxies ordered by stellar mass 
to haloes ordered by mass or circular velocity at fixed cumulative 

* E-mail: behroozi@stsci.edu 


number density provides a remarkably close match to galaxy au¬ 
tocorrelation functions, conditional stellar mass functions, satellite 
fractions, and weak lensing measurements from z = 0 to z ~ 5 (see 
also Nagai & Kravtsov 2005; Conroy et al. 2006; Watson et al. 
2015). Performing this matching at several different redshifts al¬ 
lows inferring average galaxy star formation rates and histories as a 
function of host halo mass and redshift (Conroy & Wechsler 2009; 
Firmani & Avila-Reese 2010; Leitner 2012; Bethermin et al. 2012; 
Wang et al. 2013; Moster et al. 2013; Behroozi et al. 2013b,c; 
Mutch et al. 2013; Lu et al. 2014c). This method has shown that 
average galaxy growth rates have tracked average dark matter halo 
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growth rates (multiplied by a halo-mass dependent efficiency) re¬ 
markably well over the past 12 Gyr (Behroozi et al. 2013b). 

It is less clear how closely individual galaxy growth histo¬ 
ries track individual halo growth histories (Lu et al. 2014b; Genel 
et al. 2014), especially for central galaxies (i.e., galaxies whose host 
haloes are not in orbit around any more massive halo). For exam¬ 
ple, galaxy specific star formation rates (SSFRs) at z ~ 0 show a 
clear bimodality between star-forming (SSFR > 10^" yr^^) and 
quenched galaxies (e.g., Brinchmann et al. 2004; Salim et al. 2007). 
However, the fraction of star-forming versus quenched galaxies 
falls with increasing stellar mass, whereas the fraction of host 
haloes accreting (versus losing) mass rises with increasing halo 
mass. As a result, quenching cannot be due to lack of cosmolog¬ 
ical accretion alone. Yet, it is still possible that galaxy quenching 
correlates with halo mass accretion history. For example, empirical 
models relating galaxy quenching to halo age (Hearin & Watson 
2013; Hearin et al. 2014a; Watson et al. 2015)— with older, earlier- 
forming haloes being assigned galaxies with lower star formation 
rates—have been very successful at matching quenched versus star¬ 
forming correlation functions, weak lensing, and radial profiles 
near clusters and groups. In the simplest theoretical model, a galaxy 
which uses up or expels gas faster than its host halo accretes it will 
quench due to lack of fuel (Feldmann & Mayer 2015, and refer¬ 
ences therein). Alternate scenarios could include black hole feed¬ 
back which correlates with the merger history of the halo (Silk & 
Rees 1998; Springel et al. 2005; Croton et al. 2006; Somerville 
et al. 2008), or quenching due to gravitational heating from merg¬ 
ers or infalling clumps (Cox et al. 2004; Dekel & Birnboim 2008; 
Khochfar & Ostriker 2008; Johansson et al. 2009, 2012; Birnboim 
& Dekel 2011; Moster et al. 2011). 

These quenching models have very different behaviours dur¬ 
ing major halo mergers (here, when a halo’s virial radius contains 
a smaller halo with a mass ratio of > 1 : 3). For a starvation model, 
incoming lower-mass star-forming galaxies would transfer their gas 
reservoir to the larger host, and therefore would rejuvenate star 
formation in quenched hosts. In contrast, for a gravitational heat¬ 
ing model, an incoming major merger (and/or associated accre¬ 
tion) would disrupt the flow of gas and suppress star formation. 
Finally, for a merger-fed black hole feedback model, no significant 
change would be expected until the galaxies themselves merge (as 
opposed to the smaller halo simply coming within the virial radius 
of the larger halo). These differences should be especially appar¬ 
ent for central L* galaxies (here, galaxies with 1O'®M0 < M* < 
IO'^'^Mq) at z ~ 0. Below this stellar mass range, most galaxies 
are star-forming, and above it, most galaxies are quenched (Brinch¬ 
mann et al. 2004; Salim et al. 2007), so L* galaxies are important 
probes of the quenching process. 

In this paper, we develop a galaxy pair-based selection method 
to preferentially identify haloes undergoing major mergers, and we 
examine the impact on central galaxies’ star formation rates. While 
many existing studies have found star formation enhancement in 
close pairs (e.g., Lambas et al. 2003; Nikolic et al. 2004; Alonso 
et al. 2004; Woods et al. 2006; Perez et al. 2009; Woods et al. 2010; 
Barton et al. 2007; Lin et al. 2007; Li et al. 2008; Ellison et al. 2008; 
Rogers et al. 2009; Robaina et al. 2009; Wong et al. 2011; Scudder 
et al. 2012; Xu et al. 2012; Robotham et al. 2013; Patton et al. 2013; 
Ellison et al. 2013; Scott & Kaviraj 2014), the effect size depends 
strongly on the selection process. Previous studies have typically 
focussed on pre-merging galaxies instead of merging haloes, and 
so have usually excluded the more distant pair candidates included 
in this paper (see, however, Nikolic et al. 2004; Lin et al. 2007; Li 
et al. 2008; Robaina et al. 2009; Patton et al. 2013). We also de- 



Figure 1. The stellar mass—halo mass relation at z = 0.1, from Behroozi 
et al. (2013c). A given galaxy stellar mass ratio only corresponds to the 
same host halo mass ratio in a small window near M* = lO***'^ M 0 , as seen 
by comparison to the red line with M* M),. Focussing on this mass range 
minimizes halo mass biases when selecting pairs based on stellar mass; this 
mass range is also where galaxies transition to being primarily quenched. 


sign our method to avoid past selection biases. Well-known biases 
include whether the pairs are in a cluster environment (Barton et al. 
2007), and whether the galaxies are required to be star-forming. 
Many subtler selection biases also exist. For example, the presence 
of a close companion can bias the distribution of host halo masses 
for galaxies in close pairs versus those not in close pairs; indeed, 
this has been exploited to constrain the mass of the Milky Way’s 
halo from its satellite distribution (Busha et al. 2011a; Cautun et al. 
2014). We therefore construct several mock catalogs from simula¬ 
tions for the purposes of testing for and avoiding such biases in our 
method. 

We divide the results into several sections. The observational 
selection method is described in §2, and we describe the obser¬ 
vational data sets, the mock observational catalogs, and validation 
tests in §3. We present our main findings in §4, discuss the impact 
of these results in §5 and conclude in §6. Throughout this work, 
we adopt a flat ACDM cosmology (Q.m = 0.27, flyy = 0.73,/i = 
OJ,ns = 0.95, CTg = 0.82) in close agreement with recent WMAP9 
cosmology constraints (Hinshaw et al. 2013). Stellar masses and 
star formation rates come from methods in Kauffmann et al. (2003) 
and Brinchmann et al. (2004), respectively, updated for the SDSS 
DR7, and are renormalized to a Chabrier (2003) initial mass func¬ 
tion. Halo masses are defined according to the virial spherical over¬ 
density criterion of Bryan & Norman (1998). 


2 METHOD 

Observational galaxy pair selection requires a compromise between 
simplicity, bias avoidance, and sample size. Traditionally, paired 
galaxies are selected based on being within a specified mass or lu¬ 
minosity ratio as well as within a specified projected distance and 
redshift window. Satellite galaxies in cluster environments often 
match these criteria, yet they are subject to very different physical 
conditions; many studies therefore also adopt an exclusion crite¬ 
rion, e.g., that no larger galaxy exists within a specified projected 
distance and redshift window of the galaxy pair. For simplicity and 
readers’ familiarity, we use a similar approach, but we adjust se¬ 
lection parameters so as to preferentially select major halo mergers 
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and to minimise biases as compared to previous techniques (see 
§5.5 for discussion). 

We find that a suitable selection exists for galaxy pairs where 
the larger galaxy has mass IO'^Mq < M* < 1 O**''^M 0 , and where 
the smaller galaxy’s mass is within a ratio of 0.5 dex. To simplify 
further discussion, we call the larger galaxy the host galaxy, and 
the smaller galaxy the paired galaxy. We adopt a projected radius 
cut for close pairs of 200 kpc. While larger than typical for past 
studies, this distance corresponds to the virial radius of the small¬ 
est host haloes expected to be in our sample. We also adopt a red- 
shift window of 500 km s^' for the paired galaxy, corresponding 
to the escape velocity of the largest host haloes expected to be in 
our sample. Our results are not sensitive to these specific param¬ 
eter choices, as we verify by checking many alternate choices in 
Appendix A. 

For the satellite exclusion criteria, we require that no galaxy 
more massive than the host galaxy be present within a projected 
radius of 500 kpc and a redshift window of 1000 km s^*. From 
tests with our mock catalogs (§3.2), we find that this cut retains 
77% of the central galaxies in our stellar mass range, with a purity 
of 97.4%; this is very consistent with the expected completeness 
of 75% from Liu et al. (2011). We again explore several different 
choices for these criteria in Appendix A to verify that they do not 
affect our results. 

These selection criteria are designed primarily to avoid halo 
mass biases. Dark matter haloes are roughly self-similar in their 
subhalo mass ratio distributions; however, a given ratio in dark 
matter masses for close companions can imply a very different 
ratio in their stellar masses (Fig. 1). This is why, for example, 
galaxies smaller than 10 *® M© have few massive satellites, whereas 
galaxies larger than 10 ** M 0 often have massive companions (e.g., 
Bundy et al. 2009). Selecting close pairs based on a fixed stellar 
mass ratio will therefore tend to bias the host halo masses of the 
close pair galaxies to be higher than the non-close pairs. Flowever, 
somewhat counteracting this effect, selecting close pairs within a 
fixed projected radius will probe a smaller fraction of the halo 
radius in more massive haloes. In our host stellar mass window 
( 1 O*®M 0 < Af* < 10 *®'^ M 0 ), these two biases nearly cancel each 
other out, and stellar mass ratios of 0.5 dex correspond very nearly 
to 0.5 dex ratios in halo mass—i.e., major mergers (see validation 
tests in §3.3). 

Unlike many previous studies, the host and paired galaxies 
are drawn from a stellar mass-complete sample with no require¬ 
ments on star formation activity. We also select close pairs based 
on a fixed stellar mass ratio, rather than a fixed luminosity ratio. 
Since star-forming galaxies at fixed stellar mass are brighter than 
quiescent galaxies, using a fixed luminosity ratio means that a star¬ 
forming satellite may be selected as being in a close pair, whereas 
a quiescent satellite would not be. Because of galactic conformity 
(i.e., star-forming galaxies have larger fractions of nearby star¬ 
forming galaxies than quiescent galaxies; Weinmann et al. 2006; 
Kauffmann et al. 2013; Phillips et al. 2014; Hearin et al. 2014b), 
using a fixed luminosity ratio would result in close pairs being arti¬ 
ficially more star-forming than galaxies without close companions. 

In summary, we make the following cuts for eligible host 
galaxies in our host sample: 

(i) Stellar mass between 1 O*®M 0 and 10*®'^ M 0 , and 

(ii) No more massive galaxy within 500 kpc in projected dis¬ 
tance and 1000 km s^* in redshift. 

For each of these host galaxies, we calculate the nearest paired 
galaxy, which is a galaxy that satisfies 



z 


Figure 2. Top panel: Conditional probability density for absolute r-band 
magnitude as a function of M* in the SDSS (observable volume-corrected). 
The red line shows the adopted completeness cut, which is fainter than 
96 percent of spectroscopically-targeted galaxies with mass 1O®'^-1O*®M0. 
Bottom panel: Conditional probability density for M* as a function of z in 
the SDSS. The red line shows the adopted completeness cut from the top 
panel; galaxies above the red line are taken to be volume-complete. 


□ 

□ 

SDSS gri 


Ra: 219.597 


Dec: 23.334 


119”x85” 



Figure 3. SDSS gri colour image of a BCG (119”x85”, corresponding to 
72x52 kpc at z = 0.03). The spectroscopic target selection algorithm for 
the SDSS maximises the number of non fibre-collided targets. This selects 
against BCGs, whose central location would result in fibre collisions with 
many nearby satellites. As in the above example, more spectroscopic targets 
(white squares) can be obtained in a given area if satellites are targeted 
instead of the BCG. Cases where this bias impacted host galaxy isolation 
criteria were cleaned from our sample via visual inspection. 
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(i) Stellar mass 0 to 0.5 dex less than that of the host galaxy, and 

(ii) Redshift separation within 500 km s^' of the host galaxy. 

If this paired galaxy is within 200 kpc (physical projected distance) 
of the host galaxy, we call the two galaxies a close pair. 

3 DATA AND SIMULATIONS 
3.1 Observations 

We use redshifts from the SDSS (Sloan Digital Sky Survey) Data 
Release 10 (Ahn et al. 2014), which are >90 percent complete for 
galaxies brighter than r = \1.11. In addition, we use median to¬ 
tal stellar masses and total star formation rates (SFRs) from the 
MPA-JHU value-added catalog (Kauffmann et al. 2003; Brinch- 
mann et al. 2004), updated for the imaging and spectroscopy in the 
SDSS Data Release 7 (Abazajian et al. 2009). These stellar masses 
and SFRs were calculated assuming a Kroupa (2002) initial mass 
function (IMF), which we convert to a Chabrier (2003) IMF by di¬ 
viding both by a factor 1.07. In Appendix A, we also consider fibre 
star formation rates, D„(4000) indices, the effects of BPT (Baldwin 
et al. 1981) class, and using alternate stellar masses and redshifts 
from the NYU Value-Added Galaxy Catalog (NYU-VAGC; Blan¬ 
ton et al. 2005). 

The MPA-JHU catalog covers a spectroscopic area of 8032 
deg^ and 600,480 galaxy targets with nonzero stellar masses and 
redshifts z > 0.006. We find and remove duplicate targets separated 
by a projected distance less than 2 kpc (see also Ellison et al. 2008; 
Tollerud et al. 2011), leaving 592,054 galaxies. For our tests, host 
galaxies must be further than 2 Mpc from a survey boundary or 
region of significant incompleteness; we also exclude host galax¬ 
ies with redshifts z < 0.01 to avoid Hubble flow corrections (e.g., 
Baldry et al. 2012). We exclude such boundary regions (23 percent) 
from the sample selections, although we retain them for the pur¬ 
poses of counting neighbours or testing for larger nearby galaxies; 
this leaves 479,439 galaxies over 6201.8 deg^ of sky. 

Since the SDSS is magnitude-limited, cuts are needed to con¬ 
vert to a stellar mass-complete sample. Galaxies at fixed stellar 
mass have a range of mass-to-light ratios (Fig. 2, top panel). We 
find that more than 96 percent* of galaxies in the SDSS between 
10®'^ and 10*** Mq (i.e., eligible paired galaxies) satisfy: 

M,<-0.25-1.91ogio(^^) (1) 

where is the galaxy’s Petrosian r-band absolute magnitude and 
M* is its stellar mass. 

Similarly, for galaxies at redshift z, the corresponding appar¬ 
ent magnitude limit is 

r < -0.25- 1.91og,„ -f 51og,„ (2) 

where Di{z) is the luminosity distance for our adopted cosmology. 
We note that a given host galaxy can be included in our sample only 
if a paired galaxy (§2) is detectable—i.e., only if SDSS is complete 
at the same redshift to 0.5 dex less than the host galaxy’s stellar 
mass. As making a completely volume-limited catalog would then 
unacceptably reduce sample sizes, we instead weight selected host 
galaxies inversely by the observable volume for potential paired 

' After weighting by inverse observable volume according to galaxy r-band 
absolute luminosity; this prevents bias against fainter galaxies at fixed stel¬ 
lar mass. 


galaxies, obtained by inverting Eq. 2 with r = \1.11. Galaxies in¬ 
cluded in this cut—i.e., galaxies for which r < \1.11 according to 
Eq. 2— are shown in Eig. 2, bottom panel. 

Fibre collisions can significantly reduce spectroscopic com¬ 
pleteness (Strauss et al. 2002) when multiple targetings are not 
available (~ two-thirds of the SDSS footprint). Following Patton 
et al. (2013), we reduce the observable volume by a factor 3.08 for 
galaxy pairs closer than 55” on the sky; this increases the relative 
influence of such pairs in all calculations by an equal factor. We 
have also checked pairs in the NYU-VAGC, which resolves fibre 
collisions by using the closest galaxy’s spectroscopic redshift. In 
that catalog, all close pairs within 55” are 3.23 times more numer¬ 
ous than spectroscopic close pairs within 55”, corresponding to a 
very consistent weighting factor of 3.10 after correcting for the 4% 
chance projections within 55” expected from our mock catalogs. 

The SDSS spectroscopic target selection algorithm also intro¬ 
duces an important bias near clusters. As shown in Fig. 3, targeting 
a brightest cluster galaxy (BCG) can result in many fibre collisions 
with surrounding satellites, whereas targeting the satellites instead 
would result in many more available spectra in the same region. 
Since the SDSS algorithm maximises the number of non-collided 
sources, BCGs are frequently not targeted. This affects the galaxy 
isolation criteria in §2, so we have visually examined the 1,339 host 
galaxies with paired galaxies within 500 kpc and removed 89 with 
BCGs also found within this distance. 

After applying the satellite exclusion criteria in §2, as well as 
the masking cuts and completeness corrections above, 7,303 host 
galaxies with stellar masses between 10***M0 and 10***'^ M0 re¬ 
main in the host sample; redshifts range from 0.01 to 0.057, with 
a median of z = 0.03. Of these galaxies, 439 fall in the close pairs 
sample, and of these, 36 pairs are within 55” of each other. After 
observable volume and fibre-collision weightings, 6.7% of the host 
galaxies are in a close pair. Fig. 4 shows the physically closest ten 
pairs in this sample. 

3.2 Mock Catalogs 

As a basis for mock catalogs, we use merger trees generated from 
the Bolshoi simulation (Klypin et al. 2011). Bolshoi follows 2048^ 
particles (~ 8 billion) in a comoving volume 357 Mpc on a side 
from z = 80 to z = 0 using the Art code (Kravtsov et al. 1997). The 
simulation’s particle and force resolution are 1.94 x 10** M0 and 
1.4 kpc, respectively, which correspond to a minimum resolvable 
halo mass of 10***M0. Haloes were found using the Rockstar 
halo finder (Behroozi et al. 2013d), an algorithm which determines 
particle-halo membership via a six-dimensional phase-space metric 
which is particularly suited to recovering haloes in major mergers 
and at close separations (Knebe et al. 2011, 2013a,b; Onions et al. 
2012; Pujol et al. 2014). In Appendix B, we also compare to halo 
catalogs generated using the BDM halo finder (Klypin et al. 1999; 
Riebe et al. 2013), which uses a three-dimensional density-based 
algorithm to assign particles. For both halo finders, halo masses are 
calculated using the virial overdensity criterion of Bryan & Nor¬ 
man (1998). Merger trees for both catalogs were generated using 
the Consistent Trees algorithm (Behroozi et al. 2013e), which 
compares halo catalogs across multiple timesteps to repair halo 
finder inconsistencies; the algorithm yields significantly cleaner 
mass accretion histories as compared to many other methods (Sri- 
sawat et al. 2013), especially when combined with the Rockstar 
halo finder (Avila et al. 2014). 

To generate mock catalogs, we use abundance matching, 
which assigns galaxy stellar masses to haloes with the same cu- 
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Figure 4. SDSS gri colour images of the ten closest pairs in our sample. Each image is 40x40 kpc across and is centred on the more massive (host) galaxy. A 
white “Q” denotes a quenched host galaxy (SSFR < 10“** yr“*). The separations shown here range from 5-15 kpc (median 11 kpc). 


mulative number density (Nagai & Kravtsov 2005; Conroy et al. 
2006). Several halo orderings have been explored (e.g., by mass 
or maximum circular velocity, Vmax, defined as the maximum of 

^ GM(<R\ jjjg halo’s virial radius) in Reddick et al. (2013). 

Using present-day satellite halo mass dramatically underestimates 
galaxy clustering (Conroy et al. 2006; Reddick et al. 2013), as satel¬ 
lite haloes’ dark matter is stripped much more rapidly than their 
galaxies’ stars. Instead, ordering haloes hy decreasing peak histor¬ 
ical Vmax or peak historical mass give the best matches to galaxy 
autocorrelation functions and conditional stellar mass functions. 
Specifically, Reddick et al. (2013) finds that peak historical v^ax 
gives the closest match; however, as noted in Behroozi et al. (2014), 
peak Vmax is typically set during 1:5 mergers in halo mass, and 
models using time since peak Vmax as a quenching proxy do not 
reproduce quenched galaxy distributions around clusters (Behroozi 
et al., in preparation). 

Using peak historical mass avoids these latter problems, 
but then somewhat underpredicts observed galaxy autocorrelation 
functions (Reddick et al. 2013) as many satellites continue forming 
stars even after accretion onto a larger halo (Wetzel et al. 2012). We 
can approximate this physical process with a proxy that continues 
to grow for the typical galaxy’s quenching time after a halo reaches 
its peak mass. Combining the average satellite quenching time as 
a function of stellar mass inferred in Wetzel et al. (2013) with the 
stellar mass—halo mass relation in Behroozi et al. (2013c), we ob¬ 
tain an average quenching time as a function of peak halo mass 
(Mp), which is well-approximated hy a double power-law: 


fquench(^p) — 


1.584 X 10'° yr 


10““ M, 


) “Me/ 


0.37 


(3) 


For each halo in Bolshoi, we calculate its peak historical mass, as 
well as the time at which that mass was reached (tp). If the halo is 
not currently at its peak mass (i.e., fp < tnow), we randomly select 
a continued mass accretion history from Bolshoi (from tables cal¬ 
culated in Behroozi et al. 2013a) starting at mass Mp at time tp and 
ending at tp +tquench or fnow, whichever is earlier. We call the mass 
at the end of this appended history Mq\ for haloes with tp = tnow^ 
we simply set Mq = Mp. 


Abundance matching on Mq therefore approximates the con¬ 
tinued stellar mass growth expected to happen in satellite galax¬ 
ies after their accretion, so we adopt this method for construct¬ 
ing our main mock catalogs. For the haloes under consideration, 
M M 

14.3% have jf- larger than 0.1 dex, and 3.4% have jf- larger than 
0.3 dex. For comparison, we have also constructed catalogs based 
on abundance matching with many other halo proxies (including 
peak mass and peak v^ax) to demonstrate that the catalog construc¬ 
tion method does not affect our validity tests (Appendix B). For 
the source stellar mass function, we calculate observable volume- 
corrected galaxy number counts from the same region of the SDSS 
used in §3. 1 ; this process is detailed in Appendix C. We incorporate 
a log-normal scatter of 0.2 dex in stellar mass at fixed halo mass or 
Vmax using the iterative approach in Reddick et al. (2013). Finally, 
we incorporate a Gaussian scatter of 30 km s“* in halo relative 
velocities to mimic SDSS spectroscopic redshift errors. 


3.3 Validations of the Mock Catalog and Observational Cuts 

We show several comparisons between the mock catalog (§3.2) and 
the observational galaxy sample (§3.1) in Fig. 5. We find excel¬ 
lent agreement in all cases between the mock catalog and observa¬ 
tions, including the stellar masses of the host galaxies (nontrivial 
because of the satellite exclusion criterion), the velocity separation 
between host galaxies and the nearest paired galaxy, the distribu¬ 
tion of projected distances between host galaxies and the nearest 
paired galaxy, and the large-scale environment of host galaxies, as 
measured by the number of paired galaxies within a projected dis¬ 
tance of 0.3 - 2.0 Mpc. 

Encouragingly, the large-scale environment of hosts with 
close pairs is similar to that of the entire host sample. We find 
that, while the overall satellite fraction in the host sample is small 
(2.6%), galaxies with close pairs are slightly more likely to be satel¬ 
lites (13.3%). Most of these latter cases result from scatter in the 
stellar mass—halo mass relation, which implies that the smaller of 
the two merging haloes will occasionally contain the larger of the 
two galaxies. 

The distribution of stellar masses is also similar between close 
pairs and the host sample, although our close pair selection criteria 
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Projected Distance to Nearest Paired Galaxy [Mpc] Number of Paired Galaxies with 0.3 Mpc < Rp < 2.0 Mpc 

Figure 5. Mock catalog comparisons with the SDSS. Top-left panel: host galaxy stellar mass distributions in the host sample and in close pairs; top-right 
panel: distributions of the velocity separation between host galaxies and the nearest paired galaxies (i.e., smaller galaxies with AV < 500 km s ^ and AM* < 0.5 
dex from the host); bottom-left panel: distributions of the projected distance (Rp) between host galaxies and the nearest paired galaxy; bottom-right panel: 
distributions for the number of paired galaxies with projected distances 0.3 < Rp < 2.0 Mpc, i.e., the larger-scale environment. SDSS galaxies are weighted 
by observable volume and fibre collision rate. Errors in all cases are jackknifed; mock catalog errors are smaller due to ~ 10 times larger sampled volume. 



Halo Mass [M^] 



Number of Paired Galaxies with 0.3 Mpc < Rp < 2.0 Mpc 
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0.01 0.1 1 

Projected Distance to Nearest Paired Galaxy [Mpc] 



Figure 6. Top-left panel: host halo mass distributions of close pairs compared to all galaxies in our sample, from the mock catalog. Close pairs’ host haloes 
are 20% more massive, on average, than for the host sample (see explanation in §3.4). Top-right panel: host halo mass dependence on the distance to the 
nearest paired galaxy; bottom-left panel: host halo mass dependence on environment (number of nearby paired galaxies). Bottom-right panel: halo mass 
ratios between the host galaxy and paired galaxy. Errors in the top-left and bottom-right panels are jackknifed; errors on the other two panels are bootstrapped. 
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Figure 7. Differences between total host halo mass accretion rates for close pairs and the host sample. These accretion rates encompass mass growth from all 
sources (i.e., both clumpy and unresolved “smooth” accretion). Top-left panel: median and average (stacked) mass accretion histories for the host haloes of 
galaxies in the close pairs and host samples, from the mock catalog. These show, for haloes selected at z — 0.03, how their typical mass accretion histories 
change if a close pair is also present at z = 0.03. The large bump at ^ ~ 0.2 in the close pairs sample matches the typical redshift of first accretion for incoming 
major mergers. Top-right panel: dependence of the averaged halo accretion rate on distance to the nearest paired galaxy. Separate solid lines ai‘e shown for 
averaged accretion rates over the past 100 Myr and over the period 0 < z < 0.5. Dotted lines for each timescale show the average across all paired galaxy 
separations. Middle-left panel: probability distribution of halo accretion rates averaged over 0 < z < 0.5; middle-right panel: same, averaged over past 100 
Myr. Bottom-left panel: half-mass redshift (i.e., redshift at which haloes first obtained half their z = 0 mass) distribution for close pairs and host sample; 
bottom-right panel: same, for the 4%-mass redshift. Errors on the averages in the top-right panel are bootstrapped; errors in the middle and bottom panels are 
jackknifed. 


does have a slight bias towards larger galaxies. It is possible to 
reverse this bias with an additional weighting function for close 
pairs: 

»-(«.) = 1->•8.. (ipSiS;;) M> 


This weighting is treated as a multiplicative change to galaxies’ 
observable volumes; we present both weighted and unweighted re¬ 
sults in §4. 

We note that the mocks slightly underpredict the number of 
close pairs within 200 kpc (Fig. 5, lower-left panel); host galaxies 
in the mock catalog have close pairs 6.2% of the time as compared 
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to 6.7% of the time in observations, a 9% difference. Tests with a 
higher-resolution simulation (125 Mpc on a side, 2048^ parti¬ 
cles, each with mass 2.58 x 10^ M©) show the same slight discrep¬ 
ancy, as do tests with the Millennium-II catalog from Tollerud et al. 
(2011). Several other factors could be responsible, including errors 
partitioning light in merging galaxies (see Fig. 4), slightly increased 
binding energies for galaxies as compared to their host dark matter 
haloes, scatters larger than 0.2 dex between stellar mass and halo 
mass for the host galaxies (Busha et al. 2011b), and sample vari¬ 
ance; however, most of these factors are beyond the scope of this 
paper to address. 

3.4 Host Dark Matter Halo Masses and Formation Redshifts 

We show the host halo masses corresponding to galaxies in our 
host and close pairs samples in Fig. 6 , as taken from our mock cat¬ 
alog. Close pairs have very similar host halo masses to those for 
the host sample, with only ~ 20 % larger masses (both average and 
median). This mass bias is not affected by reweighting close pairs 
according to stellar masses (Eq. 4), since there is significant scatter 
in stellar mass at fixed halo mass. Indeed, this scatter is fundamen¬ 
tally responsible for the close pairs mass bias. Since host halo mass 
correlates with the number of satellites, a larger halo mass at fixed 
galaxy mass will result in a higher likelihood of the galaxy hav¬ 
ing a close pair. Selecting galaxies with close pairs will therefore 
automatically select slightly larger haloes (Fig. 6 , top-right panel). 

Unfortunately, both stellar mass and environment (Fig. 6 , 
bottom-left panel) correlate weakly with halo mass within our sam¬ 
ple because of the highly restricted range of halo masses being con¬ 
sidered. Weighting galaxies by stellar mass or environment to fur¬ 
ther align the halo mass distributions of close pairs and the host 
sample would ruin the existing agreement between the stellar mass 
or environment distributions (Fig. 5). Weighting galaxies in close 
pairs by their velocity separations may seem attractive, but satel¬ 
lite velocities at fixed halo mass are correlated with the assembly 
time of the halo (i.e., halo age), which would introduce an addi¬ 
tional nontrivial correlation with star formation rate (Hearin et al., 
in preparation). As the cure seems worse than the problem in this 
case, we do not add any weighting to target halo masses. 

For galaxies with close pairs, the vast majority of their haloes’ 
peak mass ratios are within 1 : 3 (Fig. 6 , bottom-right panel). We 
find that for close pairs, the paired galaxy is actually within the host 
halo’s virial radius 50.5% of the time, with the remaining cases 
due to chance projections. While no selection can be perfect, this 
represents a very strong preference for major halo mergers; in the 
host sample as a whole, the major merger fraction is only 3.2%. 

We show the corresponding halo mass accretion rates (i.e., net 
matter flux into the virial radius from both clumpy and unresolved 
“smooth” sources) in Fig. 7. The host haloes of galaxies with close 
pairs experience over 100 % higher average accretion rates from 
z = 0.4 to z = 0 as compared to host haloes for the entire host galaxy 
sample (Fig. 7, top-left panel). As noted above, close pairs reside 
in somewhat larger haloes, but this would result in an expected in¬ 
crease of only 20% in their average accretion rates (Behroozi & 
Silk 2015). The enhanced accretion rates extend over two dynam¬ 
ical times before z = 0 ; this is not only due to a range of merger 
timescales, but also due to correlated structure arriving along the 
same filament as the merging halo. 

These enhanced accretion rates are only evident when a paired 
galaxy appears within the host halo’s virial radius (Fig. 7, top-right 
panel). The position of the paired galaxy within the virial radius 
does not strongly constrain infall time (Oman et al. 2013), except 


for paired galaxies very close to their host’s centres, which cannot 
have fallen in recently (Fig. 7, top-right panel). The distributions 
of host halo accretion rates are shown in the middle panels of Fig. 
7. Since halo accretion rates on 100 Myr timescales have a stan¬ 
dard deviation of 0.45 dex (Fig. 7, middle-right panel), the distri¬ 
butions for the close pairs and host sample overlap. However, the 
differences in the medians, means, and modes of the distributions 
for close pairs versus the host sample are still large: 0.3-0.4 dex, 
depending on the statistic. This large log-normal scatter in instan¬ 
taneous accretion rates also explains why average accretion rates 
are typically a factor of two higher than median accretion rates in 
the top-left panel of Fig. 7. 

Finally, we show halo assembly times in the bottom panels of 
Fig. 7. The half-mass assembly times for the haloes in the close 
pairs sample are (by construction) extremely skewed towards low 
redshifts, reflecting recent major mergers. The median half-mass 
assembly redshift for close pairs is z = 0.83, whereas it is z = 1.23 
for the host sample; 30% of close pairs have formation redshifts 
z < 0.5, whereas only 8 % of the host sample does. However, the 
longer-term halo mass accretion histories (as probed by 4%-mass 
assembly times) overlap significantly more between the close pairs 
and the host sample. 

4 RESULTS 
4.1 All Close Pairs 

We show the distribution of host galaxies’ specific star formation 
rates in the SDSS for those in close pairs and those in our full host 
sample in Fig. 8 . Galaxies with close pairs show no strong differ¬ 
ences in SSFRs, even after reweighting to eliminate a small stellar 
mass bias (Eq. 4; Fig. 8 , top-left panel). The close pairs show a 
star-forming (SSFR> 10^** yr^*) fraction of 53%, compared to 
59% for the host sample; while this is statistically very significant 
(99.6% confidence for a lower star-forming fraction), the effect size 
is modest at best. As 50% of galaxies in close pairs are undergo¬ 
ing a major halo merger (§3.4), major halo mergers can account for 
at most a 12% reduction (to 47%) in the star-forming fraction of 
central L* galaxies. The true effect is likely smaller, as the galaxies 
with close pairs also have slightly increased host halo masses and 
satellite fractions (§3.4), which would also lead to larger quenched 
fractions (Wetzel et al. 2012; Woo et al. 2013). 

Average star formation rates may be enhanced during ex¬ 
tremely close passages (i.e., projected distances < 30 kpc, or 10- 
15% of the virial radius), but no strong evidence for it is seen here 
(see, however, §4.2). Also, no evidence exists for any changes in 
the star-forming fraction of these extremely close pairs just prior to 
merging (see also Fig. 4). We have tested resolving fibre collisions 
with the redshift of the nearest galaxy to increase sample sizes, 
which further restricts the possibility of an enhanced star-forming 
fraction (Fig. A 6 in Appendix A). We have also tested expanding 
the range of host stellar masses to range from 1 O^®M 0 to lO" M 0 , 
which expands the sample size to 25,364 host galaxies, but we do 
not find any stronger evidence for enhancement of the star-forming 
fraction at close radii. 


4.2 Star-forming versus Quenched Host Galaxies 

Several previous studies have considered star-forming host galaxies 
only. For completeness, we show the radial dependence of average 
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Figure 8. Top-left panel: Specific star formation rate distribution for the host galaxies in close pairs and in the host sample, from SDSS; the weighting removes 
a small stellar mass bias for close pairs (§3.2). Top-right panel: star-forming (SSFR > 10“** yr“*) fraction for host galaxies as a function of distance from 
the host galaxy to the nearest paired galaxy (i.e., a galaxy with AF < 500 km s“* and stellar mass between 0 to 0.5 dex less than the host galaxy). Bottom-left 
and bottom-right panels: host galaxy star formation rates and specific star formation rates, respectively, as a function of distance between the host galaxy and 
the nearest paired galaxy. Median SFRs and SSFRs have been slightly offset horizontally for clarity. Errors are jackknifed for top-left panel, bootstrapped for 
all others. 


and median SFRs and SSFRs for star-forming and quenched host 
galaxies in Fig. 9. 


First, we note strong enhancements (60-100%) in SFRs and 
SSFRs for star-forming hosts with a close pair separated by < 30 
kpc. There is also a very modest average enhancement of 9.4 +^'q% 
when a close pair appears within the host halo’s virial radius (i.e., 
at about 200 kpc). This enhancement is also visible for the un¬ 
split sample (Fig. 8), although at lower significance. The quenched 
sample shows no such enhancement, perhaps by definition; how¬ 
ever, we note that the fraction of star-forming galaxies also does 
not show any enhancement (Fig. 8). This excludes chance projec¬ 
tions as a complete explanation for quenched galaxies’ behaviour. 
If quenched galaxies with close pairs were all due to chance pro¬ 
jections, and all true (3D) close pairs rejuvenated star formation in 
their hosts, then the star-forming fraction would be gradually in¬ 
creasing with decreasing pair separation—i.e., as the fraction of 
chance projections decreases—when in fact the opposite occurs 
(Fig. 8). This may imply that star-forming galaxies (but not quies¬ 
cent ones) experience a weak accretion-associated increase in SS¬ 
FRs during major mergers (see §5). We note that while quenched 
hosts appear to show a decrease in SFRs and SSFRs at close sepa¬ 
rations (< 30 kpc), this is due to sample variance, and it is not seen 
in the larger NYU-VAGC sample (Fig. A6 in Appendix A). 


5 DISCUSSION 

In §4, we found that galaxies with close pairs have only modestly 
increased quenched fractions, that close pairs do not appear to reju¬ 
venate star formation even at extremely close separations, and that 
star-forming host galaxies with close pairs show increased specific 
star formation rates (SSFRs), especially at extremely close separa¬ 
tions. We discuss how major mergers impact galaxy quenching in 
§5.1, how halo growth correlates with galaxy growth in §5.2, the 
permanence of quenching in §5.3, comparisons with semi-analytic 
models in §5.4, and comparisons with previous literature results in 
§5.5. 


5.1 Impact of Major Mergers on Galaxy Quenching 

From §4.1, central L* galaxies with a close pair within 200 kpc have 
a slightly (6%) reduced star-forming fraction, compared to central 
L* galaxies as a whole. Based on the fraction of close pairs which 
are true (3D) major mergers, major mergers could reduce the star¬ 
forming fraction by at most 12% (§4.1). The true effect is likely 
smaller, as the slightly larger host halo masses for close pairs (Fig. 
6) and slightly larger satellite fractions (§3.3) could also result in 
very modestly increased quenching. To prevent confusion, we note 
again that a major halo merger means only that the smaller halo has 
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Projected Distance to Nearest Paired Galaxy [Mpc] Projected Distance to Nearest Paired Galaxy [Mpc] 



Figure 9. Top-left panel: host galaxy star formation rates as a function of distance between star-forming host galaxies (SSFR > 10“'* yr“*) and the nearest 
paired galaxy. Top-right panel: same, for host galaxy specific star formation rates. Bottom panels: same as top panels, for quenched host galaxies. The low 
SFRs and SSFRs for extremely close pairs (< 30 kpc) near quenched hosts are not seen in larger samples (Fig. A6 in Appendix A). Errors are bootstrapped 
for all panels. 


arrived within the virial radius of the larger halo, not that the two 
haloes have become indistinguishable. 

If clump-induced gravitational heating were a major trigger 
for quenching, it is surprising that the quenched fraction changes 
so little. Close pairs are associated both with major mergers and 
enhanced accretion rates (Fig. 7, top-left panel); the peak of the 
merging activity occurs at z ~ 0.2, about one dynamical time ago. 
This is exactly the same timescale on which quenching would be 
expected to occur after suddenly shutting off cooling—for good 
reason: as both processes are driven by gravity, the characteristic 
timescales are expected to be very similar. Without a mechanism to 
delay the cooling shutoff for several dynamical timescales, gravita¬ 
tional heating is likely not an initial trigger for quenching, although 
it may help sustain pre-existing quenching. 

It also seems unlikely that major deficits in the host halo’s hot 
gas reservoir are the initial cause for quenching. Mass loss from 
infalling haloes accelerates rapidly once they pass within the virial 
radius of the host halo (e.g., Tormen et al. 1998; Kravtsov et al. 
2004; Knebe et al. 2006); hot gas can also be stripped by ram pres¬ 
sure (Kimm et al. 2011; Bahe et al. 2013). 70% of the paired galax¬ 
ies are star-forming (independent of projected distance), meaning 
that ample supplies of hot (and cold) gas should be available dur¬ 
ing major halo mergers (Popping et al. 2014; Ellison et al. 2015). 
Quenched L* galaxies are also surrounded by significant amounts 
of cold gas (Thom et al. 2012; Tumlinson et al. 2013). These results 
suggest that quenching in central galaxies may be the result of cold 


gas collapse stalling outside the galaxy (e.g., Thom et al. 2012), 
rather than a lack of gas within the halo. 


If this is the case, then tidal torques from passing galaxies may 
accelerate (and perhaps rejuvenate) star formation as the forces 
help funnel cold gas to galaxy centres (Mihos & Hernquist 1994; 
Barnes & Hernquist 1996; Hopkins et al. 2008; Patton et al. 2013). 
This hypothesis is partially consistent with the results in this study. 
The presence of an extremely close pair appears to strongly en¬ 
hance the star formation rate in star-forming galaxies (Fig. 9), but 
not to change the star-forming fraction (Figs. 4, 8, and A6). This 
may imply that the cold gas around quenched galaxies is too far 
removed from the disk to be brought in by tidal torques, or that 
the mechanism which quenches L* galaxies is strong enough to de¬ 
lay cold gas accretion until at least the physical merger of the two 
galaxies. That said, future studies using photometric stellar masses 
for extremely close pairs may be able to test for a weaker effect than 
can be constrained with our analysis. Unfortunately, using merger 
features to study this process is difficult, because morphological 
disturbances last for different amounts of time depending on the gas 
contents, morphologies, and trajectories of the progenitors (Lotz 
etal. 2011). 


© 0000 RAS, MNRAS 000, 000-000 






















Star Formation During Major Halo Mergers 11 




Figure 10. Distributions of specific star formation rates in all host galaxies and galaxies with close pairs, using identical selection as in Fig. 8 on mock catalogs 
from two semi-analytical models (Lu et al. 2012; Croton et al. 2006). Also shown are SSFR distributions for host galaxies with close pairs drawn so as to have 
identical joint probability distributions for host halo mass and stellar mass as the host sample (“Close Pairs, Mass Matched”), as well as haloes undergoing 
true (3D) major mergers (“3D Mergers, Mass Matched”). These latter two distributions allow disentangling the effects of halo mass biases in the close pairs 
sample from the effects of major mergers. 


5.2 The Relation between Halo Mass Assembly and Galaxy 
Star Formation 

We have found that, for L* galaxies, a significant recent halo merger 
event does not imply significant recent galaxy star formation. In 
addition, a lack of halo accretion does not quench star formation in 
L* and smaller satellite galaxies for several Gyr (Wetzel et al. 2013; 
Wheeler et al. 2014; Phillips et al. 2015). The galaxy’s halo mass 
and the recent assembly history (e.g., half-mass formation scale) 
therefore are not sufficient information to predict galaxy star forma¬ 
tion rates. This suggests that, at low redshifts, L* galaxy evolution 
may have become decoupled from halo evolution. In this scenario, 
the depth of the host halo’s potential well would limit the galaxy’s 
average star formation rate, but further growth of the halo would 
not correlate with growth of the galaxy.^ If true, this would imply 
that the similarity between the recent time evolutions of galaxy star 
formation rates and average halo mass accretion rates (Behroozi 
et al. 2013b) is a coincidence. At the same time, it is clear that 
stellar mass alone cannot predict quenching in central galaxies: the 
weighted close pairs sample has an identical stellar mass distribu¬ 
tion as the host sample, but a statistically significant difference in 
its quenched fraction (Fig. 8, top-left panel). ^ 

In the age-matching framework of Watson et al. (2015), star 
formation rates are correlated with halo assembly time (as mea¬ 
sured either through concentration or mass accretion history); how¬ 
ever, no previous simple measure of assembly time accounts for 
all the effects we observe. Halo concentrations (used in Watson 
et al. 2015) have the issue that major mergers have significantly 
increased concentrations during close passages (Behroozi et al. 
2014). As Watson et al. (2015) place quenched galaxies in more 
concentrated haloes, this model cannot reproduce star formation 
enhancements in very close pairs (< 30 kpc. Figs. 8 and 9). Using 
halo ages (e.g.. Fig. 7, bottom panels) would also not be able to 
reproduce strong star formation enhancements in very close pairs. 


^ This is independent of pseudo-evolution (Diemer et al. 2013), as major 
mergers contribute mass at all radial scales. 

^ While it may be concerning that the close pairs sample has a ~10% 
higher satellite fraction than the host sample (§3.3), making the isolation 
criteria stricter (for example) does not reduce the difference in quenched 
fractions (Fig. Al, Appendix A). 


However, we note that SSFR enhancements for star-forming 
hosts with close pairs are hard to interpret because of uncertainty 
in how quenched and star-forming galaxies populate haloes. While 
star-forming hosts undergoing major mergers could have up to a 
20% boost to their specific star formation rates, it is not clear 
whether this is directly attributable to the major merger, or whether 
the presence of a major merger correlates with differences in the 
longer-term mass accretion history (Fig. 7, bottom-right panel). 


5.3 The Permanence of Qnenching 

The lack of increase in the star-forming fraction for close pairs 
within 30 kpc (Figs. 8 and A6) suggests that close pairs cannot reju¬ 
venate star formation prior to merging. Nonetheless, Fig. 4 suggests 
that quenching may not be permanent. That is, star formation in an 
incoming smaller galaxy is not quenched even at extremely close 
separations; when the two galaxies merge, the star-forming regions 
in the smaller galaxy would be transferred to the larger galaxy, con¬ 
tributing to a pseudo-rejuvenation of star formation in the merger 
remnant. 

Recent stellar population analyses have shown that many 
early-type galaxies consist of a small fraction of young stars (e.g., 
Trager et al. 2000; Schiavon 2007; Zhu et al. 2010, suggesting some 
low-level star-formation has occurred recently. Imaging observa¬ 
tions in the UV/IR confirmed that star formation indeed occurs in 
many early-type galaxies (e.g., Yi et al. 2005; Kaviraj et al. 2007; 
Salim & Rich 2010. The details, such as the extent and level of star 
formation, are not yet well determined, and the responsible mech¬ 
anisms, e.g., whether or not major mergers are important, are not 
fully understood. The SAURONIFU (Integral Field Unit) observa¬ 
tions (Shapiro et al. 2010) of 48 early-type galaxies suggest that this 
low-level star formation could be either due to a large molecular 
gas content brought in by (mostly minor) mergers (Kaviraj 2014), 
or rejuvenation within the previously quenched systems. Ongoing 
IFU surveys, such as the MANGA (Bundy et al. 2015) and SAMI 
(Groom et al. 2012) surveys, will help constrain the frequency of 
these scenarios with coverage of over 10,000 nearby galaxies. 
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Table 1. Recent findings for star formation enhancement in close pairs. 


Reference 

Host Selection 

Pair Cut 

Redshifts 

Boost at < 30 kpc 

at 30-200 kpc 

Potential Confounds 

Lambas et al. (2003) 

All 2dF lOOK 

Am* < 0.75 

z^O.l 

1.2-I.4X 

N/A 

CHESi, 

Nikolic et al. (2004) 

All DRl 

Az<2 

0.03 <z< 0.1 

2-3 X 

Yes 

CHESi, 

Woods et al. (2006) 

CfA2, Zw < 15.5 

AR<2 

0.008 < z < 0.055 

Yes 

N/A 

CHS* 

Barton et al. (2007) 

2dF, Bj <-19 

Any Host 

0.010 <z< 0.088 

Yes 

N/A 

CHS* 

Lin et al. (2007)“ 

I0‘° <M*/M0 < 10“ 

Any Host 

0.1 <z< 1.1 

2-4 X 

Yes 

H 

Li et al. (2008) 

SF, DR4 

Ar< 1.4 

0.01 <z< 0.3 

1.7-2.4 X 

No 

CSo 

Ellison et al. (2008) 

SF, DR4 

AM* < 0.3 dex 

0.01 <z<0.16 

1.1-1.7X 

N/A 

CHSo 

Rogers et al. (2009) 

E,DR6 

Any Host 

0<z<0.15 

Yes 

N/A 

CHES„ 

Robaina et al. (2009) 

M* > 1O‘°M0 

AM* < 0.6 dex 

0.4<z<0.8 

2.0-2.5X 

No 

CHESi, 

Woods et al. (2010) 

All SHELS 

AR< 1.75 

0.08 < z < 0.38 

1.4-2. lx 

N/A 

CHS/, 

Wong etal. (2011) 

SF, PRIMUS, / < 22.5 

Any Host 

0.25 < z < 0.75 

1.3x 

N/A 

CHSo 

Scudder et al. (2012) 

SF, DR7 

AM* < 0.5 dex 

0.02 <z< 0.15 

1.4-2X 

N/A 

CHS„ 

Xu et al. (2012) 

SF, lO'O'* <M*/M0 < 10“ 

AM* < 0.4 dex 

0.2<z< 1 

1-1.2X 

N/A 

s„ 

Ellison et al. (2013) 

All DR4 

AM* < 0.6 dex 

0.01 ^z^O.2 

2.2-2.8X (SF) 

N/A 

CSi, 

Patton et al. (2013) 

SF, DR7 

AM* < 1 dex 

0.02 <z< 0.2 

1.5-2.2X 

Yes 

cs„ 

Robotham et al. (2013) 

GAMA, M* > 1O‘°M0 

AM* < 0.3 dex 

0.01 <z<0.089 

l-5x 

N/A 

HE 

Scott & Kaviraj (2014) 

All DR7-I-GALEX 

Any Host 

Median z ~ 0.07 

2.4-2.5X 

N/A 

CHS/, 

This Work 

DR7, 10“* < M,/M0 < 10'°'5 

AM* < 0.5 dex 

0.01 <Z< 0.057 

1.5-2X (SF) 

l.lx (SF) 



Notes. When a reference considered multiple pair criteria, the result listed is for the most major mass ratio considered. Abbreviations: DR# = SDSS (Sloan 
Digital Sky Survey) Data Release # Spectroscopic Galaxy Sample, CfA2 = Center for Astronomy Redshift Survey 2, 2dF = 2-degree Field Galaxy Redshift 
Survey, 2dF lOOK = 2dF lOOK public release, SFIELS = Smithsonian Hectospec Lensing Survey, PRIMUS = PRIsm MUlti-object Survey, GAMA = Galaxy 
And Mass Assembly, GALEX = Galaxy Evolution Explorer, SF = star-forming galaxies only, E = elliptical galaxies only. Confounds: C = galactic Conformity 
(due to luminosity delta in pair selection, luminosity threshold for sample, or star-forming paired galaxy bias), FI = Flalo mass (due to wide selection of host 
galaxy properties and/or to allowing any host galaxy to be a paired galaxy), E = Environmental effects (no isolation or environment matching for host galaxies; 
however, according to Li et al. 2008, only minor environmental effects apply to SF-only selections), S* = bias towards Star-forming galaxies (via spectroscopic 
requirements, sample luminosity threshold, and/or incomplete volume corrections), = bias against Star-forming galaxies, S^ = only Star-forming galaxies. “: 
Lin et al. (2007) uses redshift-dependent mass thresholds in the DEEP2 Galaxy Redshift Survey. 


5.4 Comparison with Semi-Analytic Model Predictions 

We have also tested mock catalogs from two semi-analytical mod¬ 
els (Croton et al. 2006; Lu et al. 2012, 2014a), which have both 
been run on the Bolshoi simulation (Fig. 10). These models have 
been tuned to match z < 0.2 stellar mass functions (Lu et al. 2014b), 
and indeed we find that statistics such as the fraction of galax¬ 
ies with close pairs (6.7% and 7.2% for the Croton and Lu semi- 
analytical models, respectively) match the SDSS result (6.7%) ex¬ 
tremely well. These models have not been tuned to match star for¬ 
mation rates, which account for the different quenched fractions 
between Figs. 8 and 10. This remains true after simulating obser¬ 
vational uncertainties in recovering SSFRs; we do so by adding 
0.2 dex of log-normal scatter to star-forming galaxies’ SSFRs and 
drawing quenched galaxies’ SSFRs from a log-normal distribu¬ 
tion with 0.35 dex scatter centered at 10^**'* yr^* (matching the 
quenched galaxy SSFR distribution in Fig. 8). 

Despite these absolute differences, the relative impact of a 
close pair on the quenched fraction in the two models remains in¬ 
teresting. Both models exhibit 20-25% larger quenched fractions 
for galaxies in close pairs as compared to the host sample (Fig. 10). 
Encouragingly, this implies that our selection approach can recover 
differences (when they exist) in star formation activity during major 
mergers. On the other hand, this implies that both models overesti¬ 
mate quenching efficiency compared to real galaxies. 

Using semi-analytic models allows disentangling the con¬ 
tributing effects (halo mass biases, increased satellite fractions, and 
major mergers) for increased quenched fractions. The effects of 
halo mass bias can be removed by creating mass-matched samples 
of close pair galaxies (i.e., all close pair galaxies with host halo 
and stellar masses within 0.1 dex of the target galaxy in the host 
sample). Similarly, the effects of major mergers can be extracted 
by considering only 3D major mergers (i.e., where the close pair is 


within the 3D physical halo radius of the host galaxy); for clarity 
of interpretation, we also mass-match this sample to the host sam¬ 
ple galaxies’ halo and stellar mass distribution. The resulting SSFR 
distributions for these mass-matched samples are shown in Fig. 10 
for both semi-analytic models. 

In both semi-analytic models, the difference between mass- 
matched and non-mass-matched close pairs samples suggest that 
half of the increase in the quenched fraction for close pairs is 
due to larger host halo masses. The remaining increase is due 
to larger fractions of “host” galaxies being satellites in the close 
pairs sample. In the Lu et al. (2012) model, quenching is imple¬ 
mented as a cooling shutoff which is a function of the host halo 
mass only, so the entire difference between the host sample and 
the mass-matched 3D major mergers sample arises from the much 
larger satellite fraction in the mass-matched major mergers sample 
(50% versus 4% in the overall host sample)."^ In the Croton et al. 
(2006) model, higher satellite fractions also yield larger quenched 
fractions for the mass-matched major mergers sample. However, 
in this model, galaxies are quenched more indirectly, using black 
hole feedback. Since cooling is not shut off, the presence of a ma¬ 
jor merger increases the amount of gas available to cool onto the 
galaxy, increasing SSFRs in star-forming host galaxies with true 
(3D) major halo mergers by a factor of 2 compared to star-forming 
galaxies in the host sample. In the original close pairs sample, the 
Croton et al. (2006) model predicts SSFRs enhanced by 34% for 
star-forming galaxies compared to the host sample, which is sig¬ 
nificantly larger than our findings in the SDSS (enhancements of 
~ 10%; §4.2). 

This satellite fraction is much larger than in the abundance-matched mock 
catalogs, which suggest that the satellite fraction in the close pairs sample 
should be 13.3% (§3.3), and that the satellite fraction in the mass-matched 
3D major mergers sample should be 19.0%. 
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5.5 Comparison with Previous Work 

Table 1 summarises several recent literature results on star for¬ 
mation enhancements in close pairs. The majority of previous re¬ 
sults preferentially include star-forming galaxies in their samples. 
There are good reasons for doing so, especially for improved statis¬ 
tics (due to star-forming galaxies being brighter at fixed mass). 
While conformity and halo mass biases may he present in these 
past works, they appear to have a modest impact on measured en¬ 
hancements for extremely close pairs. Compared to our results, Li 
et al. (2008) found ~ 20% larger enhancements using a luminosity 
delta for extremely close pairs (resulting in a conformity bias; see 
§2) while using similar stellar masses and star formation rates from 
earlier MPA-JHU catalogs. That said, biases on the 20% level can 
be important when studying star formation enhancements at larger 
separations. 

Among these past results, there is general agreement that ex¬ 
tremely close pairs (separated by < 30 kpc) show star formation 
enhanced by 50-150%, although the reported values range from 
20% to 400%. Less agreement exists about whether enhancements 
persist at larger radii, but this may be due to methodological limita¬ 
tions. Papers which report no enhancement (Li et al. 2008; Robaina 
et al. 2009) use photometric identification of pairs; however, within 
150 kpc, the ratio of true satellites to background contaminants is 
~ 15% (Liu et al. 2011). Hence, the weak enhancements (5-20%) 
reported in the papers that use spectroscopic pair selection (Nikolic 
et al. 2004; Lin et al. 2007; Patton et al. 2013) would be reduced to 
percent-level effects in the photometric selections. While not con¬ 
sidered in this paper, past results suggest that star formation en¬ 
hancements are lower for minor mergers (Woods et al. 2006, 2010; 
Scudder et al. 2012) and for more massive galaxies (Li et al. 2008; 
Robotham et al. 2013). 

Robotham et al. (2013) is the only paper in Table 1 that also 
creates a stellar mass-limited local pair sample including both star¬ 
forming and non-star-forming host galaxies. Unfortunately, due to 
limited statistics, their constraints on star formation enhancements 
in extremely close pairs are weak (0-400%). With modestly better 
statistics, our results constrain this to be between 0% and 100%; 
with the still better statistics in Fig. A6, this becomes 50-100%. 
Robotham et al. (2013) report halo mass (FoF group-based) differ¬ 
ences between their close pairs and control samples, and also do not 
isolate their host galaxies; however, both effects would be expected 
to reduce any enhancements, rather than increase them. 


6 CONCLUSIONS 

We have presented a method for selecting central L* SDSS galaxies 
whose host haloes are preferentially undergoing major mergers. 
Tests with mock catalogs (§3.3) suggest that our selection can 
identify host galaxies with near-identical host halo masses, stellar 
masses, and environments, but with average halo mass accretion 
rates higher by 0.3 dex over the past 5 Gyr (§3.4). Additionally, 
50% of galaxies selected in this way are undergoing major halo 
mergers, as compared to 3% of isolated L* galaxies (§3.4). 

Our findings include: 

(i) The subsample with 50% major mergers has a 6% lower star¬ 
forming fraction than the whole isolated host galaxy sample, imply¬ 
ing at most a 12% effect in a pure major mergers sample. (§4.1). 

(ii) This latter finding limits how gravitational heating or gas 


reservoir transfers in mergers can affect central galaxy star forma¬ 
tion rates (§5.2). 

(iii) Consistent with previous research, star-forming host galax¬ 
ies show 70% larger SSFRs when an extremely close pair is present 
(<30 kpc), but only ~ 10% larger SSFRs when the paired galaxy 
is between 30-200 kpc in projected separation (§4.2). 

(iv) Extremely close pairs (<30 kpc) do not appear to rejuvenate 
star-formation for quenched host galaxies (§4.1). 

(v) Quenching for central L* galaxies does not depend on their 
stellar mass alone (§5.2). 

(vi) Previous halo age-based or concentration-based methods 
for matching galaxy star formation rates to haloes have a difficult 
time reproducing all SFR enhancements found for SDSS galaxies 
(§5.2). 

(vii) Current semi-analytical models over-predict the impact of 
major mergers on galaxy quenching (§5.4). 
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APPENDIX A: VARIATION IN SELECTION 
PARAMETERS 

We have tested varying the isolation criteria for host galaxies, the 
definition of a paired galaxy, as well as the separation distance con¬ 
sidered “close.” In Figs. Al and A2, we show the impact on the 
specific star formation rate distributions as well as the overall star¬ 
forming fraction of host galaxies. We find no impact large enough 
to affect our main conclusions; e.g., that major halo mergers do not 
result in significant changes to star formation activity. 

We have also explored two spectroscopic-only indicators of 
galaxy formation. These have not been included in the main dis¬ 
cussion because the SDSS fibre size (3” diameter) only covers a 
small fraction of the host galaxy, which may not be representa¬ 
tive of the overall star-formation activity (Salim et al. 2007). The 
fibre-only SSFRs (Brinchmann et al. 2004) shown in Fig. A3 do 
not provide any different picture than total galaxy SSFRs (Figs. Al 
and A2). The 4000A break strength (D„(4000)) is more interesting. 
As with galaxy total SSFRs, a larger fraction of host galaxies with 
close pairs seem to have older stellar populations (D„(4000)> 1.7). 


However, star-forming host galaxies in the close pairs sample ap¬ 
pear to have slightly younger ages than star-forming galaxies in the 
host sample. This difference is statistically significant, but small 
(median D„(4000) of 1.359^QQgg for host galaxies with close pairs, 
versus 1.397 )')qqq 2 for the host sample). While this could indicate a 
younger stellar population, changes of this magnitude are also pos¬ 
sible with dust and metallicity differences (Heman-Caballero et al. 
2013). Future IFU spectroscopy of close pairs (e.g., Bundy et al. 
2015) or deeper spectroscopy at higher redshifts may reveal larger 
differences in galaxy discs. 

Many previous works have excluded galaxies classified as 
AGNs (based on the BPT diagram; Baldwin et al. 1981) due to 
the difficulty of estimating star formation rates from emission lines. 
The Brinchmann et al. (2004) SFRs avoid this issue by basing AGN 
and composite fibre SFRs on D„(4000), using the D„(4000)-SSFR 
distribution for non-AGN-contaminated galaxies. While this could 
introduce small biases for AGN and composite galaxies, the low 
redshifts of our sample (z < 0.06) mean that SDSS fibre sizes (3”) 
cover a minority of the galaxy light. As a result, the majority of 
the SFR estimate is based on photometry, which is calculated in the 
same way for all galaxies. Nonetheless, for completeness, we show 
galaxy SSFR distributions for the host and close pairs samples in 
Fig. A5 excluding AGNs, LINERs, and composite host galaxies. 

Finally, we have tested using stellar masses and redshifts 
from the NYU Value-Added Galaxy Catalog (NYU-VAGC; Blan¬ 
ton et al. 2005). This catalog’s main advantage is that redshifts for 
fibre-collided galaxies are taken from the nearest available galaxy. 
This more than triples the available statistics for close pairs sepa¬ 
rated by <55”. However, as in Fig. 8, there remains no evidence 
for extremely close pairs having a larger star-forming fraction (Fig. 
A6, top panels). Similarly, boosted specific star formation rates are 
seen in star-forming host galaxies in close pairs, whereas no boost 
is seen for quenched galaxies, regardless of pair separation (Fig. 
A6, lower panels). 


APPENDIX B: ALTERNATE MOCK CATALOG 
CONSTRUCTION METHODS 

As noted in Reddick et al. (2013), many different ways exist to 
abundance match galaxies to haloes; additionally, many different 
halo finders exist (see Knebe et al. 2011, 2013b, for comparisons 
and reviews). We have therefore explored three alternate ways of 
assigning galaxies to haloes in mock catalogs. These include abun¬ 
dance matching on Mp (peak historical mass) and Vjnax.peak (peak 
historical v^ax) with the Rockstar halo finder, as well as abun¬ 
dance matching on Mp with the BDM halo finder. 

As shown in Fig. Bl, the choices of abundance matching 
proxy and halo finder both affect quantitative details for galaxy 
host halo masses (or Vmax) and accretion history. Several qualitative 
facts remain, however. Regardless of the catalog, our selection cri¬ 
teria for close pairs does not significantly bias present-day host halo 
properties relative to those for the host sample. The major merger 
fractions for the close pairs samples are 47%, 52%, and 44% for 
the Rockstar Mp, Rockstar v^ax.peak- and BDM Mp catalogs, 
respectively. Additionally, regardless of the catalog, host haloes of 
galaxies in close pairs have significantly larger recent total accre¬ 
tion rates as compared to the host sample. In combination, these 
suggest that the selection criteria we have chosen are a robust way 
to preferentially identify galaxies whose host haloes are undergoing 
major mergers and have had more recent formation times without 
imposing a strong bias on the host halo mass or Vmax- 
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Figure Al. Specific star formation rate distributions for all host galaxies and those with close pairs when selection criteria are varied. For reference, the default 
selection criteria for host galaxies are: no larger galaxy within 1000 km s“^ in redshift (“Isolation Av”) or 0.5 Mpc in projected distance (“Isolation Distance”); 
for a smaller galaxy to be called a paired galaxy, it must be within 500 km s“* in redshift (“Paired Galaxy Av”) and 0.5 dex in stellar mass (“Paired Galaxy 
AM*”). In all panels, close pairs are defined as a paired galaxy within 200 kpc in projected distance from a host galaxy. For alternate definitions of close pair 
distances, see Fig. A2. 
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Projected Distance to Nearest Paired Galaxy fMpc) 



Projected Distance to Nearest Paired Galaxy [Mpc] 
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Figure A2. Star-forming fraction for host galaxies as a function of distance to the nearest paired galaxy when selection criteria are varied, as in Fig. Al. 
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Figure A3. Left panel: distribution of fibre SSFRs for all host galaxies and those in close pairs. Right panel: fibre SSFRs as a function of distance to the nearest 
paired galaxy. The diameter of an SDSS fibre is 3”, corresponding to 0.43 kpc at z = 0.01 and 2.4 kpc at z = 0.057. For our stellar mass range, typical host 
galaxy half-light radii are 1-2.5 kpc (Kravtsov 2013). 




D^(4000) Projected Distance to Nearest Paired Galaxy [Mpc] 

Figure A4. Left panel: distribution of Drt(4000)—i.e., 4000A break strength—for all host galaxies and those in close pairs (Balogh et al. 1999 definition). A 
weak 4000A break (D„(4000)< 1.6) signifies a young stellar population, and a strong one (D„(4000)> 1.7) signifies an older stellar population. Right panel: 
D„(4000) as a function of distance to the nearest paired galaxy. The diameter of an SDSS hbre is 3”, con'esponding to 0.43 kpc at z = 0.01 and 2.4 kpc at 
z = 0.057. For our stellar mass range, typical host galaxy half-light radii are 1—2.5 kpc (Kravtsov 2013). 




Figure A5. Left panel: distribution of SSFRs for all host galaxies and those in close pairs, excluding hosts classified as AGN or LINERs according to the BPT 
diagram. Right panel: same, except also excluding composite host galaxies. 
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Figure A6. Top-left panel: the star-forming fraction of host galaxies as a function of the distance to the nearest paired galaxy, using stellar masses and redshifts 
from the NYU-VAGC (Blanton et al. 2005). In this catalog, galaxies missing spectroscopic redshifts are assigned the redshift of the nearest neighbour galaxy, 
which results in a ~ 3 times larger sample size for close galaxy pairs within 55”. This panel is analogous to Fig. 8, top-right panel. Top-right panel: SSFRs 
for host galaxies (still from Brinchmann et al. 2004) as a function of distance to the nearest paired galaxy, analogous to the bottom-right panel of Fig. 8, 
using NYU-VAGC stellar masses and redshifts. Bottom-left panel: SSFRs for star-forming host galaxies as a function of distance to the nearest paired galaxy, 
analogous to the top-right panel of Fig. 9, using NYU-VAGC stellar masses and redshifts. Bottom-right panel: SSFRs for quenched host galaxies as a function 
of distance to the nearest paired galaxy, analogous to the bottom-right panel of Fig. 9, using NYU-VAGC stellar masses and redshifts. 


APPENDIX C: STELLAR MASS EUNCTION 
CALCULATION 

The local volume (z < 0.07) is underdense compared with the 
nearby (0.07 < z < 0.2) universe (see Baldry et al. 2008; Keenan 
et al. 2013, and references therein). When abundance-matching to 
a simulation at the universe’s typical density, it is necessary to 
correct for this effect; otherwise the stellar mass-halo mass rela¬ 
tion and satellite fractions will be underestimated. The local cos¬ 
mological underdensity results in a relatively uniform reduction in 
galaxy counts at all masses (Baldry et al. 2008), so massive galaxy 
counts can be used to trace the underdensity as a function of red- 
shift. Number densities for M* > 10 *'Mq galaxies as a function 
of redshift in our SDSS sample are therefore shown in Fig. Cl. We 
model the local underdensity as linearly dependent on redshift out 
to z = 0.0644, where we assume that it reaches the cosmological 
mean (see fit in Fig. Cl). The corresponding correction factor for 
galaxy number counts is then: 


In the absence of detailed completeness information as a function 
of redshift and Kauffmann et al. (2003) stellar mass, we have also 
scaled all number counts to account for the average SDSS spec¬ 
troscopic completeness fraction of 92%. The resulting stellar mass 
function (both before and after corrections) is shown in Fig. Cl, 
with comparisons to previous literature results (Li & White 2009; 
Moustakas et al. 2013). 
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^maxpeak ^ Average Halo Accretion Rate for 0<z<0.5 [M^ yr 

Figure Bl. Left panels: Peak halo mass (or peak v^ax^ as appropriate) for haloes hosting galaxies in the host sample as well as the close pairs sample for three 
alternate mock catalogs. Right panels: averaged halo accretion rates for 0 < z < 0.5 for the same alternate mock catalogs. 
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Figure Cl. Left panel: number density of > 10^^ M© galaxies in the SDSS as a function of redshift. Horizontal error bars show bin widths; data point 
centres are located at the median bin redshift. Vertical error bars show Poisson uncertainties. Right panel: corrected and uncorrected stellar mass functions 
compared to previous literature results. 
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